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1. Introduction 

HERA is the only electron(positron)-proton collider in the world. It has 
been running since 1992, with a major upgrade of the accelerator during 
2001-2002. The multipurpose HI and ZEUS detectors have been collect- 
ing data steadily during these years. Pre-upgrade integrated luminosities of 
~ 100 pb^ 1 per experiment have already allowed the most accurate deter- 
mination of the proton structure up to date. Both collaborations have also 
made important contributions in making tests and precision measurements 
of QCD and the electroweak sector of the Standard Model, search for new 
particles and new interactions, study of jet production and jet substruc- 
ture, study of the structure of the photon, heavy flavour production, and 
investigation of diffraction and the structure of the pomeron. 

This report deals with the tests of QCD on hadronic final states, concen- 
trating mainly on the results from jet production. Multijet QCD production 
will be the main background for searches of Higgs, supersymmetric particles 
and new physics in general. Also, the luminosity of the colliding particles 
(partons from the protons) at LHC will be governed by QCD and so precise 
determinations of the protons parton distribution functions (PDFs) and the 
strong coupling constant, a s , will be crucial for any cross section measure- 
ment. In this sense, ep collisions at HERA provide a suitable enviroment 
to do precision studies of QCD in hadronic-induced reactions (as opposed 
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to e + e~ annihilations at LEP) in a cleaner set-up than in pp collisions at 
Tevatron. 

The main sources of jets at HERA are neutral current (NC) deep inelas- 
tic scattering (DIS), with Q 2 ^> Aq CD , where Q 2 is the square of the mo- 
mentum transfer, and photoproduction (PHP), with Q 2 ~ 0. Up to leading 
order (LO) in a s , jet production in NC DIS proceeds via the quark-parton 
model (Vq — > q, where V = 7* or Z°), boson-gluon fusion (Vg — » qq) and 
QCD-Compton (Vq — > qg) processes. The jet production cross section is 
given in perturbative QCD (pQCD) by the convolution of the proton PDFs 
and the partonic cross section. Only two kinematic variables are needed to 
characterise these processes; they can be taken as Q 2 and Bjorken x. The 
inelasticity variable y, which gives the energy transfer in the proton rest 
frame, is also useful. In PHP, a quasi-real photon emitted by the electron 
beam interacts with a parton from the proton to produce two jets in the 
final state. In LO QCD, there are two processes which contribute to the jet 
photoproduction cross section: the resolved process, in which the photon 
interacts through its partonic content, and the direct process, in which the 
photon interacts as a point-like particle. The jet production cross section is 
given by the convolution of the flux of photons in the electron, the parton 
densities in the proton and photon and the partonic cross section. Thus, 
measurements of jet cross sections in NC DIS and PHP can be used to per- 
form tests of pQCD, help to constrain the gluon density in the proton and 
determine a s . 



ZEUS 




(a) (b) (c) 

Fig. 1. (a) Inclusive-jet cross section in NC DIS as a function of Q 2 pQ; (b) inclusive- 
jet cross section in NC DIS as a function of Ej' B in different regions of Q 2 0; (c) 
energy-scale dependence of a s pQ. 
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Fig. 2. (a) Thrust distribution in different Q 2 regions [3]; (b) ao and a s (Mz) values 
extracted from event-shape distributions 4 . 
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Fig. 3. Forward-jet cross section as a function of Borken x 5 compared to (a) LO 
and NLO calculations; (b) Casacade MC predictions; and (c) CDM and Rapgap 
MC predictions. 



2. Jet search in ep collisions 

Jets in hadronic collisions are not as easily identified as in e + e~ annihi- 
lations since the initial-state colliding partons carry only a fraction of the 
energy of the parent hadrons and are accompanied by several soft hadrons 
not correlated with the hard interaction. The spectator partons give rise 
to the so-called remnant jets (hadronisation of spectator partons), underly- 
ing event (soft interactions between the spectator partons) and multiparton 
interactions (hard interactions between the spectator partons). 
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Fig. 4. Azimuthal correlation as a function of x for different Q 2 regions |S] compared 
to the CDM and Casacade MC predictions. 
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Fig. 5. (a) x° hs distribution for data and resolved, direct and resolved+direct MC 
events [7]; (b) Dijet cross section as a function of | cos#*| in PHP for ir° bs ^ 0.75 [7|. 



The collisions in hadronic-type reactions do not occur in the centre-of- 
mass frame; the initial-state partonic system is boosted along the beam axis 
by an amount which is different for every event. To treat on equal footing all 
possible final-state hadronic systems, variables invariant under longitudinal 
boosts, such as transverse energy (Et), differences in pseudorapidity (tj) 
and azimuth ((f>), are best suited to reconstruct jets in this type of reactions. 
Also, the use of the transverse energy helps to identify the signal for a hard 
scattering and to disentangle the products of the hard interaction from the 
beam remnants. These features apply to both NC DIS and PHP processes at 
HERA, though in NC DIS there is a further complication due to the fact that 
the exchanged boson is virtual and carries pt- This effect is compensated 
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Fig. 6. (a) Inclusive-jet cross section as a function of (b) a s as a function 

of By* 0; (c) inclusive-jet cross sections as a function of if et [§]. 
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Fig. 7. Three-jet cross sections in PHP for M 3 ' > 25 GeV as functions of (a) x° bs 
and (b) y For comparison, the predictions of Pythia and Herwig with and 
without MPIs are included. 



for by selecting a frame in which the virtual boson and the proton collide 
head-on. This is the so-called Breit frame. 

At HERA, the Ut cluster algorithm in the longitudinal inclusive mode 
has been proven to be the best algorithm to reconstruct jets in these hadronic- 
type reactions. This algorithm has been used by both the HI and ZEUS 
collaborations since many years for making presicion tests of pQCD. The 
advantages of this algorithm can be summarisedas follows: (i) it allows a 
transparent translation of the experimental set-up to the theoretical calcu- 
lations (it avoids the ambiguities related to overlapping and merging of jets 
present in cone algorithms); (ii) pQCD calculations using this algorithm 
are finite to all orders (as opposed to the cone algorithm, for which the 
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Fig. 8. (a) Inclusive-jet cross section in NC DIS as a function of £y B in different 
regions of Q 2 £Q; (b) Dijet cross section in PHP as a function of E 1 ^" in different 
regions of rf ct for x° bs > 0.75 0. Both measurements are compared to NLO 
calculations based on ZEUS-jets proton PDFs 13 . 
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Fig. 9. (a) Valence, sea and gluon distributions from the ZEUS-jets proton PDF 
parametrisations (b) uncertainties of the ZEUS-jets proton PDF parametri- 
sations [T5| . 



calculations diverge for orders higher than the next-to- leading) ; and (m) 
it presents the smallest hadronisation corrections compared to other algo- 
rithms. This algorithm has been used to reconstruct the jets in data and 
theoretical predictions for the measurements presented in the next sections. 
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Fig. 10. (a,b) Summary of a s (Mz) measurements at HERA compared with the 
world and HERA averages JTj; (b) combined measurements of a s (E^ t ) from 
HERA El- 



3. QCD studies in neutral current deep inelastic scattering 

3.1. Inclusive-jet cross sections and a s 

Inclusive-jet cross sections have been measured in NC DIS at large Q 2 
(Q 2 > 125 GeV 2 [J and Q 2 > 150 GeV 2 as functions of Q 2 , the jet 
transverse energy in the Breit frame (J?y' B ) and for Ej, B in different regions 

of Q 2 (see Fig. ^) . The measured cross sections as functions of Q 2 and Ej, B 
display a steep fall-off of several orders of magnitude within the measured 
range; the data show a harder spectrum in E^ B as Q 2 increases. The 
measurements have been compared to next-to-leading-order (NLO) QCD 
calculations. The predictions give a good description of the data over a 
wide range of Q 2 and £f B . This proves the validity of the description 
of the dynamics of inclusive-jet production by pQCD at (D(aa 2 ). The LO 
prediction for the jet cross section in the Breit frame is directly proportional 
to ot s . Values of a s (Mz) have been extracted from the measurements: 
a s (M z ) = 0.1212 ±0.0017 (stat.) ±g;8oif ( ex P-) -oS ( th -) (Q 2 > 500 GeV 2 ), 
a s (M z ) = 0.1186 ± 0.0030 (exp.) ± 0.0051 (th.) (Q 2 > 150 GeV 2 ), 
The small experimental (2.9% and 2.5%, respectively) and theoretical 
(2.3% and 4.3%, respectively) uncertainties make these determinations the 
most precise at HERA and comparable to those values extracted from more 



8 



plhc2006proc printed on February 7, 2008 



inclusive measurements. The energy-scale dependence of a s has been tested 
from the measurements of the inclusive-jet cross section as a function of 
E^ B . The results, shown in Fig. ^c), are in very good agreement with the 

running of a s as predicted by pQCD over a wide range in Ejf-g. 

3.2. Event shapes and the hadronisation process 

The hadronic final-state in NC DIS has also been used to study the 
hadronisation process, which is a non-perturbative effect. Recent develop- 
ments in the model of power corrections permit the understanding of these 
effects from first principles. The model of Dokshitzer et al. is characterised 
by an effective universal coupling, 6lq. Event-shape variables, like thrust, 
broadening, C parameter and jet mass (inspired by e + e~ measurements), 
are well suited to test the universality of this effective coupling. In this 
type of analysis, the data are compared to a model prediction which con- 
sists of a combination of NLO QCD calculations and the expectations of 
the power-correction model: F = Fp er turbative + -Fpower correction, where F 
is an event-shape mean or distribution. In the case of distributions, the 
pertubative prediction is supplemented by matched resummed calculations. 
Figure [2[a) shows, as an example, the measured thrust distribution in dif- 
ferent regions of Q 2 j3j together with the predictions which have been fitted 
to the data to extract cxq and a s (Mz). It is possible to obtain a good de- 
scription of the event-shape observables in hadronic-induced reactions in the 
regions of phase space where the predictions of the power-correction model 
are valid. The results of the fits to extract «o and a s (Mz) are shown, for all 
event-shape observables measured, in Fig. |2fb) [1]. The extracted a s (Mz) 
values are consistent for all observables and with the world average (shaded 
band). It is possible to extract a universal non-perturbative parameter, 
a.Q = 0.5 ± 10%, from all the measured observables. This supports the 
concept of power corrections as an appropiate alternative approach for the 
description of hadronisation effects. 

3.3. Parton evolution at low x and unintegrated PDFs 

One of the main channels of Higgs production at LHC is expected to 
be gg — > H. Predictions for this process need information on the parton 
evolution at low x and on the unintegrated proton PDFs. Forward-jet data 
at HERA are ideally suited to study these effects. 

At high scales (Q, E^), NLO calculations using the DGLAP evolution 
equations provide a good description of the data (see Section 3). Thus, mea- 
surements at HERA have provided accurate determinations of the proton 
PDFs and precise determinations of a s . The DGLAP evolution is equiva- 
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lent to the exchange of a parton cascade in which the exchanged partons 
are strongly ordered in virtuality kx- This approximation works well at 
high scales, as mentioned above, but it is expected to break down at low x 
since it includes only resummation of leading logarithms in Q 2 and neglects 
contributions from logl/x, which are important when logQ 2 <C logl/x. 
Approaches to parton evolution which deal with low x comprise: (i) BFKL 
evolution, which includes resummation of logl/x terms and exhibits no &t 
ordering; (ii) the CCFM equations, which incorporate angular-ordered par- 
ton emission and are equivalent to BFKL for x — > and to DGLAP at 
large x; (Hi) the introduction of a second fc^-ordered parton cascade on 
the photon side a la DGLAP, which is implemented by assigning a partonic 
structure to the virtual photon, mimicks higher-order QCD effects at low 
x. HERA is an ideal testbed for these theoretical approaches. These tests 
are realised by restricting the jet data to large jet pseudorapidity values 
(forward direction, proton side), Xj et = E^ /E p 3> XBj (to suppress quark- 
parton- model type events) and Q rs Ej, (to restrict evolution in Q 2 ). 

FigureOJshows the measured forward-jet data |S] as a function of Bjorken 
x. The data increases as xsj decreases. The data are compared to differ- 
ent predictions. The NLO pQCD calculations, based on DGLAP evolution 
equations, fail to describe the data at low x. The predictions of the Cas- 
cade Monte Carlo (MC), which is an implementation of the CCFM equa- 
tions, give an improved description of the data at low x, but the description 
at high x worsens. The predictions of the CDM model, which contain no-/cr 
ordering, and the resolved-photon MC give a better description of the data, 
except at very low x. 

Insight into low-x dynamics can also be gained by studying the azimuthal 
separation (A^>) between the two hardest jets. An excess of events at small 
A(f> would signal a deviation from DGLAP evolution. The fraction S of dijet 
events with < 120 has been measured as a function of Bjorken x in 
different regions ofQ 2 . The data are shown in Fig.|I]compared with different 
predictions. The predictions of the CDM model give a good description of 
the data at low x and low Q 2 , but are below the data for high Q 2 . The 
predictions of the Cascade MC using different sets of unintegrated PDFs 
are also shown: the predictions display sensitivity to the unintegrated gluon 
distributions. Therefore, these measurements can be used to constrain the 
unintegrated PDFs. 

4. QCD studies in photoproduction 

4-1- Test of color dynamics in dijet PHP 

Measurements of jet cross sections in PHP allow tests of color dynam- 
ics. At HERA, quark and gluon exchange can be studied in dijet PHP by 
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separating resolved and direct processes using the fraction of the photon en- 
ergy invested in the production of the dijet system, which is measured via 
the observable x° bs = ^2 E^e -11 ^ 1 /2yE e . For direct events, x° bs ~ 1 and 
for resolved events, x° bs < 1. Resolved processes are dominated by gluon 
exchange (like dijets in pp collisions) and the angular behaviour of the cross 
sections is oc (1 — |cos#*|)~ 2 for 9* — > (Rutherford scattering), where 9* 
is the scattering angle in the dijet centre-of-mass frame. Direct processes 
proceed via quark exchange (similar to prompt-photon processes in pp col- 
lisions) and the angular behaviour of the cross section is oc (1 — | cos#*|) -1 
for 9* — > 0. Therefore, the cos9* distribution reflects the underlying parton 
dynamics since it has sensitivity to the spin of the exchanged particle in 
two-body processes. 

The data sample for x° bs > (<)0.75 is expected to be dominated by 
direct (resolved) processes, as demonstrated by the x° bs distribution for data 
and resolved, direct and resolved -l-direct MC events shown in Fig. Efa) [7j. 
These samples were used to compute the cross section as a function of 
| cos 9*\ (see Fig. Efb)) 0. The cross sections increase as |cos6**| increases, 
but the cross section for x° bs < 0.75 displays a much more rapid increase 
for | cos (9*| — > 1 than the cross section for x° bs > 0.75. The LO and NLO 
QCD calculations, which include resolved (direct) processes with an angular 
behaviour oc (1 — | cos 0*[) -2 ^ 1 ) give a very good description of the data. 
This demonstrates that the underlying parton dynamics in PHP is well 
understood. 



4-2. Inclusive-jet cross sections, a s and the photon PDFs 

Measurements of jet cross sections in PHP allow tests of pQCD. Inclusive- 
jet cross sections have been measured in PHP at large E^ (E^ > 17 GeV jH] 
and Ei^ 1 > 21 GeV (U]) as functions of and rf et (see Figs. EJa) and (c), 
respectively). The measured cross section as a function of Ej, displays a 
steep fall-off of several orders of magnitude within the measured range. The 
measurements have been compared to NLO QCD calculations. The predic- 
tions give a good description of the data. This proves the validity of the 
description of the dynamics of inclusive-jet production by pQCD at 0{aa 2 s ). 
The QCD prediction for the jet cross section is directly proportional to a s . 
A value of a s {Mz) has been extracted from the measurement of da/dE^: 
a s (M z ) = 0.1224 ± 0.0001 (stat.) t°;°° 22 (exp.) ±g;ggg (th.). 

The small experimental (1.8%) and theoretical (4.4%) uncertainties make 
this determination also one of the most precise at HERA. This value and 
the ones quoted in Section 3.1 are in agreement with the world |10j and 
HERA averages. The energy-scale dependence of a s has been tested 
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from the measurements of the inclusive-jet cross section as a function of 
. The results, shown in Fig. [BJb), are in very good agreement with the 
running of a s as predicted by pQCD over a wide range in 

Measurements of jet cross sections in PHP also allow tests of the photon 
PDFs. The structure of the photon can be investigated at HERA by mea- 
suring jet cross sections most sensitive to the photon PDFs, namely dcr/drf ct 
(see Fig. Etc)) 9 or do~/dx° bs , and comparing the measurements to predic- 
tions based on different parametrisations of the PDFs. The measurements 
can then be used to discriminate among different parametrisations or used 
in a global fit to constrain them, as it was recently done for the proton 
PDFs (see Section 5). 

4-3. Multijet cross sections and multiparton interactions 

Measurements of jet cross sections in PHP allow tests of hard multijet 
production and multiparton interactions (MPIs), which are expected to be 
copious at LHC energies. Multijet production is directly sensitive to higher 
orders since e.g. the three-jet cross section is proportional to a 2 s at lowest 
order. Such measurements allow tests of the parton showers in the MC 
models. Their sensitivity to MPIs allows tests and tuning of the models. 
Figure [7| shows the three-jet cross section for M 3 J > 25 GeV as functions of 
x° bs and y ^] • The data are compared to the predictions of the P ythia 
and Herwig MC models which describe the shape of the y distribution but 
fail to describe the shape of the measured x° bs distribution. The inclusion 
of MPIs in Pythia, tuned to generic collider data, results in a failure to 
describe both data distributions. In the case of Herwig, the MPI model 
was tuned to the x° bs distribution and gives a good description of the data, 
but the description of y gets spoiled. Therefore, these very precise data are 
an ideal ground for tuning and testing the parton-shower and MPI models 
in a clean hadronic-induced enviroment. 

5. Jet cross sections and the proton PDFs 

Measurements of jet cross sections in NC DIS and PHP are directly sensi- 
tive to the proton PDFs and provide a useful constrain on the gluon density, 
especially at mid- to high-x values, which are most relevant at LHC ener- 
gies. This is achieved in the case of PHP by restricting the measurements 
to regions of phase space less sensitive to the photon PDFs, i.e. x° bs > 0.75. 
Recently, very precise measurements of inclusive-jet cross sections in NC 
DIS ciS cl function of E^ B in different regions of Q 2 (see Fig. Eta)) Hj and 

dijet cross sections in PHP as a function of Ej, in different if et regions for 
x° bs > 0.75 (see Fig. Efb)) QJ] have b een incorporated in a rigorous way in 
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a QCD fit to determine the proton PDFs [T^j. The valence, sea and gluon 
distributions as a function of x for Q 2 = 10 GeV 2 are shown in Fig. Eta). 
The Incorporation of jet cross sections in the fit resulted in an improve- 
ment on the determination of the gluon density in the proton: as shown in 
Fig-EJb), the uncertainty in the gluon density decreases for mid- to high-x 
by up to a factor of two. 

6. Conclusions 

HERA has become an unique QCD-testing machine, very useful for un- 
derstanding multijet production in clean hadronic-induced reactions. Con- 
siderable progress in understanding and reducing uncertainties has led to 
very precise determinations of a s (see Fig. HOf a)). The dominant uncer- 
tainty in these measurements comes from the thoretical side thus, improved 
QCD calculations are needed for better accuracy. The running of a s has 
been observed from HERA jet data alone (see Fig. llOf b)), New HERA I re- 
sults are still coming; new determinations of a s (Mz) are shown in Fig. llOf c) 
and more than 400 pb _1 of integrated luminosity of e^p data with polarised 
lepton beams are being analysed. Therefore, there is a wealth of hadronic- 
type-reaction data already available to test new MC models, MC@NLO and 
next-to-next-to-leading-order pQCD calculations if they are made available 
for ep collisions. 
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